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ABSTRACT 

We use high-resolution dissipationless simulations of the concordance flat ACDM model to make 
predictions for the galaxy-mass correlations and compare them to the recent SDSS weak lensing mea- 
surements of Sheldon et al. (2004). The simulations resolve both isolated galaxy-size host halos and 
satellite halos (subhalos) within them. We use a simple scheme based on the matching of the circular 
velocity function of halos to the galaxy luminosity function and on using the observed density-color 
correlation of the SDSS galaxies to assign luminosities and colors to the halos. This allows us to closely 
match the selection criteria used to define observational samples. The simulations reproduce the ob- 
served galaxy-mass correlation function and the observed dependence of its shape and amplitude on 
luminosity and color, if a reasonable amount of scatter between galaxy luminosity and circular velocity 
is assumed. We find that the luminosity dependence of the correlation function is primarily determined 
by the changing relative contribution of central and satellite galaxies at different luminosities. The color 
dependence of the galaxy-mass correlations reflects the difference in the typical environments of blue and 
red galaxies. We compare the cross-bias, b x = b/r, measured in simulations and observations and find a 
good agreement at all probed scales. We show that the galaxy-mass correlation coefficient, r, is close to 
unity on scales > lh^ 1 Mpc. This indicates that the cross bias measured in weak lensing observations 
should measure the actual bias b of galaxy clustering on these scales. In agreement with previous studies, 
we find that the aperture mass-to-light ratio is independent of galaxy color in the range of luminosities 
probed by observational samples. The aperture mass scales approximately linearly with luminosity at 
L r > 10 10 ft.~ 2 L©, while at lower luminosities the scaling is shallower: A/as oc L® 5 . We show that 
most of the luminous galaxies (M r < —21) are near the centers of their halos and their galaxy-mass 
correlation function at r < lOO/i -1 kpc can therefore be interpreted as the average dark matter density 
profile of these galaxies. Finally, we find that for galaxies in a given narrow luminosity range, there is 
a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation 
between mass and luminosity derived from the weak lensing analyses should be interpreted with caution. 



Subject headings: cosmology: theory 
of universe 

1. INTRODUCTION 

Understanding the processes that shape the clustering of 
dark matter and galaxies is one of the main goals of obser- 
vational cosmology. Modern large redshift surveys, such 
as the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS, Colless et al. 2001) and the Sloan Digital Sky Survey 
(SDSS, York et al. 2000), allow measurements of galaxy 
and galaxy-mass correlations and of their dependence on 
galaxy properties and environment with unprecedented ac- 
curacy. 

Concurrently, cosmological A-body simulations have de- 
veloped into a powerful tool for calculating the gravita- 
tional clustering of collisionless dark matter in hierarchi- 
cal cosmologies with well-specified initial conditions. The 
main obstacle in direct comparisons between models and 
data is understanding the dependence of theoretical pre- 
dictions on both the relatively straightforward physics of 
gravitational clustering and the more complex physics of 
galaxy formation. The former determines the distribution 
of dark matter, while the latter affects the relationship 
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between the distribution of galaxies and mass, the "bias". 
The complexity of processes operating during galaxy for- 
mation on a very wide range of scales makes it difficult to 
include them directly in simulations, although efforts in 
this direction are ongoing (e.g., Katz et al. 1999; White 
et al. 2001; Yoshikawa et al. 2001; Pearce et al. 2001; 
Berlind et al. 2003; Weinberg et al. 2004). 

Given the difficulty of large-scale galaxy formation sim- 
ulations, two simpler approaches have recently been pur- 
sued to make theoretical predictions for galaxy proper- 
ties. In the first, the galaxies are identified with dark 
matter halos and subhalos in dissipationless cosmological 
simulations (e.g., Colin et al. 1999; Kravtsov & Klypin 
1999; Neyrinck et al. 2003; Kravtsov et al. 2004a). In the 
second hybrid approach, collisionless A-body simulations 
are combined with a semi-analytic treatment of galaxy 
formation (e.g., Kauffmann et al. 1997; Governato et al. 
1998; Kauffmann et al. 1999a,b; Kolatt et al. 1999; Benson 
et al. 2000a,b; Somerville et al. 2001; Wechsler et al. 2001; 
Berlind et al. 2003). 

Galaxy-galaxy lensing is a relatively new but important 
observational probe of the relation between galaxies and 
dark matter. It directly measures the galaxy-mass cross- 
correlation function around galaxies of different types and 
environments (see, e.g., Brainerd et al. 1996; Fischer et al. 
2000; Hoekstra et al. 2001, 2004; McKay et al. 2001; Shel- 
don et al. 2001, 2004; Smith et al. 2001; Wilson et al. 
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2001). 

Several theoretical studies of the galaxy-mass correla- 
tions in cosmological simulations have been carried out 
in the last three years (Guzik & Seljak 2001, 2002; Yang 
et al. 2003; Weinberg et al. 2004). Still, there appear to 
be several discrepancies between the various model results 
and observations, which are yet to be properly understood. 
Guzik & Seljak (2001) used intermediate-resolution dark 
matter only simulations combined with a semi-analytic 
model to identify galaxies and make predictions for the 
galaxy-mass correlation function. They compared these 
results to the first weak lensing detection from SDSS (Fis- 
cher et al. 2000), and found the amplitude of the galaxy- 
mass correlation in their model to be systematically higher. 
They attributed this discrepancy to the differences be- 
tween the luminosity function used in their semi-analytic 
prescription and the observed luminosity function. Yang 
et al. (2003), using the same theoretical model (intermediate- 
resolution dark matter simulations combined with the Kauff- 
mann et al. 1999a semi-analytic galaxy formation pre- 
scription), compared to more extensive observational re- 
sults from McKay et al. (2001). They also found that 
the galaxy-mass correlation function in the simulations 
was systematically higher than the observational measure- 
ments, by about a factor of two. Consequently, they found 
mass-to-light ratios about a factor of two higher than those 
observed. The study of Weinberg et al. (2004) investigated 
several of the same statistics, based on galaxies identified 
in SPH simulations. They found that for systems with 
small baryonic masses, their dark matter-to-baryon ratios 
agreed with the mass-to-light ratios, derived from weak 
lensing data, presented by McKay et al. (2001). For large 
masses, however, the baryonic masses of their simulated 
galaxies were high by a factor of 1.5-2 compared to what 
is required for agreement with the SDSS data. In each 
of these studies, discrepancies with the data may have 
been due to discrepancies of the model luminosity func- 
tions with that observed. Guzik & Seljak (2002) used a 
more phenomenological approach, and developed a formal- 
ism in the context of the halo model to fit the McKay et al. 
(2001) data. They extracted information on, e.g., the de- 
pendence of the virial mass-to-light ratio on luminosity, 
the typical mass for galaxies, and the fraction of galaxies 
in groups and clusters. 

The most recent observational study of galaxy-galaxy 
lensing by the SDSS collaboration (Sheldon et al. 2004, 
hereafter S04), has significantly improved the accuracy 
of galaxy correlation measurements, and may shed light 
on many of the previous uncertainties. The spectroscopic 
sample of 127,001 lensing galaxies was combined with more 
than 9 x 10 6 source galaxies with photometric redshifts, 
which allowed a detailed study of the relation between 
mass and light for several luminosity bands, morpholog- 
ical types, and on scales from 20 h^ 1 kpc to 10 h^ 1 Mpc 
(physical). In addition to verifying the findings of previous 
studies, S04 identified new features in the data, including 
a scale-dependent luminosity (cross) bias. 

Besides being a general test of the Cold Dark Matter 
paradigm, comparisons with cosmological simulations can 
be used to gain a deeper insight into the interpretation 
of these new observational results. Furthermore, they can 
shed light on to what extent weak lensing results can be 



Table 1 
Simulation parameters 



Name a$ L box N p m p /i peak 

/r'Mpc h- 1 M Q h- 1 kpc 



ACDM 6 oa 


1.0 


60 


256 3 


1.07 x 10 9 


1.9 


ACDM 6 ob 


0.9 


60 


256 3 


1.07 x 10 9 


1.9 


ACDM 6 oc 


0.9 


60 


256 3 


1.07 x 10 9 


1.9 


ACDMgoa 


0.75 


80 


512 3 


3.16 x 10 8 


1.2 


ACDMgob 


0.9 


80 


512 3 


3.16 x 10 8 


1.2 


ACDM120 


0.9 


120 


512 3 


1.07 x 10 9 


1.8 



used to estimate halo masses, or to learn about the aver- 
age dark matter halo density profiles. In this paper we use 
high-resolution dissipationless simulations of the concor- 
dance ACDM model to study the galaxy-mass correlation 
function with specific emphasis on comparing to the ob- 
servational measurements of S04. Results of several recent 
studies suggest that gravitational dynamics is the domi- 
nant mechanism shaping galaxy clustering, at least in the 
simple case of galaxies selected above a luminosity or mass 
threshold (Kravtsov & Klypin 1999; Kravtsov et al. 2004a; 
Zcntner et al. 2004). Kravtsov et al. (2004a), using dark- 
matter only simulations which resolve galactic mass sub- 
halos, matched galaxies of a given luminosity to a popula- 
tion of halos and subhalos of a given circular velocity and 
the same number density, and found excellent agreement 
with the galaxy-galaxy correlation functions measured in 
the SDSS. This study extends that approach to investigate 
the galaxy-mass correlations, and makes detailed compar- 
isons to the new Sheldon et al. (2004) results. 

The paper is organized as follows. In § 2 and § 3 we 
describe the simulations and the halo identification algo- 
rithm. The halo samples used in our analysis are described 
in § 4. The main results, including a detailed comparison 
with the most recent SDSS measurements, arc presented 
in § 7. In § 8 and § 9 we discuss and summarize our results 
and conclusions. 

2. SIMULATIONS 

For each simulation in this study we assume the concor- 
dance flat ACDM model: S\> = 1 - ^A = 0.3, h = 0.7, 
where f2 and Oa are the present-day matter and vac- 
uum densities, and h is the dimensionless Hubble constant 
defined as Ho = lOO/i km s _1 Mpc -1 . This model is con- 
sistent with recent observational constraints (e.g., Spergel 
et al. 2003; Tegmark et al. 2004). To study the effects of 
the power spectrum normalization, the effects of box size, 
and cosmic variance, we consider a set of simulations listed 
in Table 1. The simulations are labeled according to their 
box size. 

Three different simulations of 60/i _1 Mpc s=a 85.71 Mpc 
box with 256 3 ~ 1.67 x 10 7 particles are used. One of 
these three simulations is normalized to erg = 1.0, where 
erg is the rms fluctuation in spheres of 8/i _1 Mpc comov- 
ing radius. This simulation was used previously to study 
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the halo clustering and bias by Kravtsov & Klypin (1999) 
and Colin et al. (1999); the reader is referred to these 
papers for further details. The other two 60/i _1 Mpc sim- 
ulations have a lower normalization of ag = 0.9, and differ 
only in the realization of the Gaussian initial conditions. 
The second subset of simulations followed the evolution of 
512 3 « 1.34 x 10 8 particles in a Mpc « 114.29 Mpc 

box. These simulations have different power spectrum nor- 
malizations of a s — 0.75 and ag = 0.9, but were started 
from the same realization of the initial conditions. We also 
use a simulation of a 120 h~ x Mpc box with 512 3 particles 
and ag — 0.9. 

The simulations were run using the Adaptive Refine- 
ment Tree iV-body code (ART; Kravtsov et al. 1997; Kravtsov 
1999). The ART code reaches high force resolution by 
refining all high-density regions with an automated re- 
finement algorithm. The criterion for refinement is the 
mass of particles per cell. In the ACDMeoa the code starts 
with zeroth level uniform 512 3 grid and refines an individ- 
ual cell only if the mass exceeds n t h — 5 particles, inde- 
pendent of the refinement level. In terms of overdensity, 
this means that all regions with overdensity higher than 
<5 = fith 2 31 /n, where n is the average number density of 
particles in the cube, are refined to the refinement level 
L. Thus, for the ACDM 6 o a simulation , n is 1/8. The 
peak formal dynamic range reached by the code in this 
simulation is 32, 768, which corresponds to a peak formal 
resolution (the smallest grid cell) of ft pea k = 1.83ft- -1 kpc; 
the actual force resolution is w 2ft pca k = 3.7/i _1 kpc (see 
Kravtsov et al. 1997). In the higher-resolution ACDM 80 
simulations the initial grid is 256 3 and the refinement cri- 
terion is level- and time-dependent. At the early stages of 
evolution (o < 0.65) the thresholds are set to 2, 3, and 4 
particle masses for the zeroth, first, and second and higher 
levels, respectively. At low redshifts, a > 0.65, the thresh- 
olds for these refinement levels are set to 6, 5, and 5 par- 
ticle masses. The lower thresholds at high redshifts are 
set to ensure that collapse of small-mass halos is followed 
with higher resolution. The maximum achieved level of 
refinement is L max = 8, which corresponds to a comov- 
ing cell size of 1.22/i _1 kpc. As a function of redshift the 
maximum level of refinement is equal to L max = 6 for 
5 < z < 7, L m ax = 7 for 1 < z < 5, L max > 8 for z < 1. 
The peak formal resolution is /i poa k < \.2h~ x kpc (physi- 
cal). The refinement criteria for the ACDM120 simulation 
are similar to those of the ACDMso runs, except that ini- 
tially the entire volume is resolved with a 1024 3 grid. The 
parameters of the simulations are summarized in Table 1. 

3. HALO IDENTIFICATION 

3.1. The algorithm 

A variant of the Bound Density Maxima halo finding 
algorithm (Klypin et al. 1999) is used to identify halos 
and the subhalos within them. The details of the algo- 
rithm and parameters used in the halo finder can be found 
in Kravtsov et al. (2004a). The main steps of the algo- 
rithm are identification of local density peaks (potential 
halo centers) and analysis of the density distribution and 
velocities of the surrounding particles to test whether a 
given peak corresponds to a gravitationally bound clump. 
More specifically, we construct density, circular velocity, 
and velocity dispersion profiles around each center and it- 



eratively remove unbound particles using the procedure 
outlined in Klypin et al. (1999). We then construct final 
profiles using only bound particles and use these profiles to 
calculate properties of halos, such as the circular velocity 
profile V c i rc (r) = y/GM(< r)jr and compute the max- 
imum circular velocity V max . Our completeness limit is 
approximately 50 particles. This limit corresponds to the 
mass below which cumulative mass and velocity functions 
start to flatten significantly (Kravtsov et al. 2004a). In 
the following analysis, only halos with masses M > 50m p 
are considered. 

3.2. Halo Classification: host halos, satellites, and 
central galaxies 

In this study we distinguish between host halos with cen- 
ters that do not lie within any larger virialized system and 
subhalos located within the virial radii of larger systems. 
Below the terms satellites, subhalos, and substructure are 
used interchangeably. To classify the halos, we calculate 
the formal boundary of each object as the radius corre- 
sponding to an enclosed overdensity of 180 with respect 
to the mean density of the universe. Halos whose centers 
are located within the boundary of a larger mass halo are 
classified as subhalos or satellites. 

The halos that are not classified as satellites are iden- 
tified as host halos. Note that the center of a host halo 
is not considered to be a subhalo. Thus, host halos may 
or may not contain any subhalos with circular velocity 
above the threshold of a given sample. The host centers, 
however, are included in clustering statistics because we 
assume that each host harbors a central galaxy at its cen- 
ter. Therefore, the total sample of galactic halos contains 
central and satellite galaxies. The former have the posi- 
tions and maximum circular velocities of their host halos, 
while the latter have the positions and maximum circular 
velocities of subhalos. 

In the observed universe, the analogy is simple. In a 
cluster, for example, the brightest central galaxy that typ- 
ically resides near the center would be associated with the 
cluster host halo in our terminology. All other galaxies 
within the virial radius of the cluster would be considered 
"satellites" associated with subhalos. 

4. THE GALAXY SAMPLES 

In this study, our galaxy sample is created by assign- 
ing realistic SDSS luminosities and colors to dark matter 
halos. To construct mock galaxy catalogs for comparison 
with observations, one must define selection criteria for 
particular halo properties to mimic the selection function 
of the observational sample as closely as possible. Halo 
mass is often used to define halo catalogs: e.g., a catalog 
can be constructed by selecting all halos in a given mass 
range. However, the mass and radius are poorly defined 
for the satellite halos due to tidal stripping — which alters 
a halo's mass and physical extent (see Klypin et al. 1999). 
Therefore, we use the maximum circular velocity, Vmax, as 
a proxy for the halo mass. For isolated halos, Vmax and 
the halo's virial mass are directly related. For subhalos, 
V max will experience secular decrease but at a relatively 
slow rate (Kravtsov et al. 2004b). 

To mimic the observational selection function, r-band 
luminosities are assigned to the halos as follows. We match 
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the cumulative velocity function n(> V max ) of the halos to 
the SDSS observed r-band cumulative luminosity function 
(Blanton et al. 2003c). Note that n(> V max ) includes both 
isolated host halos and subhalos. We use the r-band data 
since it is in that band that the SDSS spectra were se- 
lected, in which the luminosity function is more reliably 
measured observationally, and on which most SDSS anal- 
yses focus. The mean rcdshift of SDSS galaxies in the 
main sample is z — 0.1, and this is also quite close to the 
weighted mean for the lens galaxies used by S04. We thus 
start by matching the M r (z = 0.1) luminosity function 
to the halos identified at z — 0.1, obtaining the average 
V max — M r relation by matching n(> V max ) to n{< M T ). 
This was the same method used to assign luminosities to 
halos in Kravtsov et al. (2004a), in which galaxy cluster- 
ing properties were reproduced remarkably well. To match 
what was done with the S04 magnitudes, we A-correct the 
z = 0.1 magnitudes to z = (using the code kcorrect 
v3_2, described in Blanton et al. 2003a). 

Note that the magnitudes in S01 have not been corrected 
for evolution effects, and we use the Blanton et al. (2003c) 
luminosity function appropriate for z = 0.1. The bright- 
est galaxies in the S01 sample are typically at a somewhat 
higher rcdshift. Therefore, if we had included evolution 
effects the galaxy-mass correlation of our brightest sample 
would be lower than is shown below, since the brightest 
objects would then be too bright to be included in the 
corresponding S01 sample. However, the effect is quite 
small in the magnitude range we consider: smaller than 
the uncertainty due to cosmic variance and luminosity as- 
signment for galaxies. The evolution correction for the 
brightest galaxies is less than 0.1 mags, and would result 
in a less than 10% correction to the AS amplitude. 

Observations indicate that there is substantial scatter 
in the relation around the mean Tully-Fisher (TF) and 
Faber-Jackson correlations for spiral and elliptical galax- 
ies. The magnitude of the scatter varies for different galaxy 
samples, wavebands, slopes of the correlation, the correc- 
tions applied, etc. Willick et al. (1997) found that the 
TF residuals are consistent with a Gaussian distribution 
with the standard deviation equal to the rms scatter after 
various corrections to the absolute magnitudes are applied. 
The rms values for this scatter in magnitude at a given ve- 
locity in red and infrared wavelengths range from <~ 0.2 to 
<~ 0.9 magnitudes (e.g., Aaronson & Mould 1983; Willick 
et al. 1997; Verheijen 2001). 

Motivated by these findings, we introduce scatter in the 
relation between V max and M r , assuming that at a given 
Vmax there is a Gaussian distribution of possible M r val- 
ues. The luminosities are assigned to halos using an it- 
erative procedure that preserves the observed luminosity 
function. We start by deriving the mean M r — V maK rela- 
tion by matching the cumulative velocity function in the 
simulation with the observed luminosity function. Given 
our assumption of a Gaussian scatter in M r with fixed 
standard deviation um t or, equivalently, in natural log lu- 
minosity with a fixed standard deviation <j\ n l , the scatter 
in velocity at fixed luminosity will be given by <j\ n v max = 
dlogVmax I dlogL x <7\ n l , where the V max — L relation is 
taken to be locally well approximated by a power law. As 
we cannot change the observed galaxy luminosities and 
luminosity function, to assign a luminosity accounting for 




Fig. 1. — The r-band absolute magnitude (the magnitudes arc 
M — 5 log h) vs. the maximum circular velocity of halos. The solid 
curve is obtained by matching the cumulative velocity function n(> 
V m ax) to the SDSS luminosity function n(< M r ), both at z = 0.1. 
The dashed line denotes the mean Vmax as a function of M r , after 
scatter between luminosity and Vmax is introduced as described in 
the text. For comparison, the dotted lines show a power law L r oc 
KSax' with a = 3.5 and a = 1. 



scatter we choose to perturb the values of halo V max us- 
ing the expression above and using the mean V max — M r 
relation. Note that the maximum circular velocities are 
perturbed only during the procedure of assigning lumi- 
nosities to halos and the actual V ma x of an object does 
not change. In each iteration, we thus perturb the ve- 
locities by a Gaussian in natural log with this standard 
deviation, construct a new circular velocity function, and 
match this to the observed luminosity function to assign 
new luminosities to halos. The above procedure gives a 
new mean M r — V~ max relation which is used to assign scat- 
ter for the next iteration. The iterations are repeated until 
convergence. 

In the analyses presented below we assume a standard 
deviation of 1.5 magnitudes for the M r distribution at 
fixed Vn iax . This magnitude of scatter is needed to match 
the amplitude of the galaxy-mass correlation for the bright 
samples with the observed values. Although this rms is 
larger than the values derived from the TF analyses, we 
note that the latter are derived for the luminosities cor- 
rected for the effects of extinction, inclination, etc. We, 
on the other hand, are applying scatter to the raw lu- 
minosities uncorrected for any such effects. For example, 
the correction due to internal extinction alone can be as 
high as one magnitude (e.g., Verheijen 2001). We make 
the simplifying assumption that the scatter is just a func- 
tion of V m ax, and do not model it as a function of type 
for a given V max . Including all galaxy types in a sam- 
ple is also likely to increase the scatter (Kannappan et al. 
2002; Courteau et al. 2003). The magnitude of scatter we 
assume is therefore perfectly consistent with the observa- 



MODELING GALAXY-MASS CORRELATIONS IN SIMULATIONS 



5 



tions. Given that the scatter does affect the amplitude of 
the galaxy-galaxy and galaxy-mass correlation for bright 
galaxies, it should not be neglected. However, since the 
fiducial value we assume is still uncertain, we note where 
our results may be sensitive to this assumption. 

The relation between the r-band absolute magnitude 
and the maximum halo circular velocity obtained with and 
without accounting for scatter is shown in Figure 1. The 
scatter has a significant effect. Without the scatter, the 
relation has a running slope. Including the scatter trans- 
forms the relation into almost a power law L r oc V^ K . 
The effect of scatter is similar to the Malmquist bias: be- 
cause the abundance of galaxies increases with decreasing 
luminosity, more low luminosity objects scatter into sam- 
ple at a fixed mass than high luminosity objects. As we 
show below, the effect of scatter on the amplitude of AS 
is strong at high luminosities and is small at L < L*. 

The limited size of the simulation box puts an upper 
limit on the luminosities of galaxies that can be studied 
reliably, and the completeness limit of our halo catalogs 
imposes a lower limit on the luminosities of objects that 
can be considered. Thus, we will not present a compar- 
ison of our results to the S04 sample as a whole, since 
our sample lacks both the faintest of their objects (—18 < 
M r < —17), and, most importantly, their brightest bin 
(-24.0 < M r < -22.2). To obtain reliable statistics for 
objects that bright in simulations, a considerably larger 
box size is needed. Instead, wc try to mimic the interme- 
diate luminosity samples of S04 as closely as possible. 

To investigate trends in other bands and with color, 
magnitudes are assigned to subhalos in the remaining four 
SDSS bands (u, g, i, and z) using the observed relation 
between local galaxy density and color (similar to the pro- 
cedure described in Wechsler 2004, Wechsler et al. 2004). 
We first select a volume-limited sample of SDSS galaxies 
from the first public data release (DR1, Abazajian et al. 
2003). We use the CMU-Pitt Value Added Catalog to 
compile a subsample of these galaxies which is not sensi- 
tive to edge effects, and to get a local density measurement 
for each galaxy: the projected distance to the tenth near- 
est galaxy neighbor, brighter than M r — 51og/i = —19.7 
and within cz = 1000 km/s. We measure the identical 
quantity for our mock galaxies. For each mock galaxy, we 
then choose a real SDSS galaxy which has a similar r-band 
luminosity and nearest neighbor distance, and assign the 
colors of this galaxy to the mock galaxy (bins are chosen to 
contain ~ 50 galaxies) . Assigning colors in this way allows 
us to test whether the color- and band- dependent trends 
in the data can be reproduced with dark matter halos, un- 
der the simple assumption that the dominant variables for 
determining galaxy colors are its luminosity (or the circu- 
lar velocity of its subhalo) and the local environment (as 
is indicated by Blanton et al. 2003b). 

5. ERROR ESTIMATES 

In what follows, all distances, both three-dimensional 
and projected, are comoving. We study the mass-mass 
and galaxy-mass two-point correlation functions at dis- 
tances from 20/i _1 kpc up to approximately 15/i _1 Mpc, 
except for the galaxy-galaxy correlation function, which 
cannot be measured down to 20h~ 1 kpc in our runs. The 
simulation data is divided into 20 equal logarithmic bins in 



1 

io- 1 

10 5 
10 4 
10 3 
10 2 
10 

1 




-Smith et al. (2003)- 

EuJ I I i 



TTTJ 1 1 Mill 

cosmic variance 



n — i i i i i in 1 — i i i i 1 1 1 



r fi 



5* 



CM 



: I I hh 

uJ i I i 



13 



from top to bottom: 

^mm ' ^hm' ^hh 



- ACDM B 



= 0.9 



ACDM ao <7 B =0.75 



io- 



1 

r (h -1 Mpc) 



10 



Fig. 2. — Top left panel: The effect of the box size. Top right 
panel: The effect of cosmic variance. Bottom panel: the effect of 
power spectrum normalization. In each panel we plot, from top 
to bottom, the mass-mass (£ mm ), the halo-mass (fj, m ), and the 
halo-halo (£hh) two-point correlation functions. For all boxes, halos 
(including subhalos) of the same number density are used to calcu- 
late the correlation functions. For clarity purposes, the halo— mass 
and the mass-mass functions are scaled up by a factor of 10 and 
100, respectively. The jackknife resampling error bars are plotted. 
The parameters of the simulations are summarized in Table 1. All 
correlation functions are presented at z = 0.1. 



r. As an estimate of the error in the correlation function 
in each bin we use the maximum of the Poisson error and 
the jackknife resampling error. The latter is an estimate of 
the cosmic variance and is calculated by dividing the sim- 
ulation box in eight sub-volumes. The jackknife variance 
(Tj K of the quantity \ ' IS then 



N — 



^JK = 



N 



1 N 



Xif 



(1) 



Here, N — 8 is the number of sub-volumes, \% is the mean 
of x when the i-th sub-volume is excluded, and \i 1S the 
mean of Xi- Typically, the jackknife resampling error is 
larger than the Poisson error at all scales. 

6. COSMIC VARIANCE AND THE EFFECTS OF BOX SIZE 

AND <7 8 

In this section, results for the different simulations are 
compared to understand the effects of the box size, of dif- 
ferent power spectrum normalizations (as), and of cosmic 
variance on the two-point correlation functions. The re- 
sults are shown in Figure 2. In each panel the mass-mass, 
halo-mass, and halo-halo correlations are plotted from top 
to bottom, respectively Each is measured for a volume- 
limited sample that consists of objects of the same number 
density for all boxes. The halo mass and mass-mass cor- 
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Fig. 3. — Left panel: AE as a function of projected separation R as measured by Sheldon et al. (2004) (points) and as calculated 
from the simulations (lines). The results arc shown for two r-band luminosity bins: —22.2 < M r < —21.7 (triangles / dashed line) and 
—21.7 < M r < — 17. 0(— 18.0) (squares and solid line, respectively). The effect of its can be seen for the fainter sample for which we plot 
the results for the erg = 0.75 (dotted line) and erg = 0.9 (solid line) runs. Right panel: The same as the left panel but showing the (3D) 
galaxy-mass correlation function. The error bars plotted are 1 a jackknife resampling errors. For clarity, error bars in the simulation results 
are plotted only for scales at which there appears to be some discrepancy between simulations and observations. For illustration purposes, 
the results for the most luminous bin arc shifted upward by a factor of 10 in both panels. Also shown on the right panel with thin lines are 
the results when no scatter between the luminosity and V m ax is introduced. The effect of scatter on AS is similar. 



relation functions are displaced by a factor of 10 and 100 
upward, respectively, for clarity 

The effect of simulation box size can be seen in the up- 
per left panel of the figure. Each correlation function is 
shown here for three different box sizes, all with exactly the 
same cosmological model and power spectrum normaliza- 
tion. Previous studies based on the mass-mass correlation 
function have found that the size of the simulation box af- 
fects the amplitude on large scales (e.g., Itoh et al. 1992; 
Colombi et al. 1994; Colin et al. 1999). According to these 
studies, a small box leads to an underestimate of the cor- 
relation function on scales larger than <~ 10% of the box 
size. The finite size of the box sets an upper limit on the 
longest wavelength of fluctuations that can be present in 
the simulation. At a given separation, however, the two- 
point correlation function, takes (weighted) contributions 
from all modes. This may be the case with ACDMgo c 
and ACDMgob. From our comparison it becomes clear 
that cosmic variance can significantly alter all these ex- 
pected trends. Furthermore, while in the linear regime 
each mode evolves independently, in the nonlinear regime 
mode coupling occurs and power leaks from large to small 
scales. We check the effects of non-linearity by using the 
model developed by Smith et al. (2003) for the calcula- 
tion of non-linear mass power spectra. The prediction of 
their model is presented as the long dashed-dottcd line — 
hardly discernible in the figure, since the agreement of our 
simulations with the model is very good. 

The upper right panel demonstrates the amount of cos- 
mic variance for the box sizes we use. Simulations ACDMgofc 
and ACDM 6 o c are identical in every respect, except that 



different random realizations of the initial conditions were 
used. The comparisons show that the jackknife error may 
underestimate the cosmic variance somewhat. The effect 
of varying the value of a% is shown in the lower panel for the 
two 80 h~ x Mpc boxes. Note that these two simulations 
used the same random realization of the initial conditions. 
Since the same structures form in both 80ft. -1 Mpc simula- 
tions, all differences between the two simulations are thus 
due to the difference in cr 8 . 

Keeping in mind these points, in what follows we focus 
on two of our largest volume simulations, ACDMsoh and 
ACDMi2o- The latter provides more robust results at high 
luminosities (M r < — 21) due to better statistics than the 
smaller volume simulations. Based on the presented box 
comparisons, the jackknife error bars should be indicative 
of how different our results can be if another simulation 
were used. 

7. RESULTS 

7.1. Weak tensing observables 

We now compare our simulation results to those recently 
measured by S04 using the SDSS. Using the shapes of the 
background galaxies, they estimated the tangential shear 
7t azimuthally averaged over thin projected radial annuli 
from the lens galaxy. The average tangential shear is given 
in terms of the quantities related to the lens mass surface 
density (Kaiser 1984; Miralda-Escude 1991). More specif- 
ically, 

_ 7t = AE/Ec-it, AE = E(< R) - E(i?) (2) 
where E(< R) is the mean surface density within the pro- 
jected radius R, E(i?) is the azimuthally averaged sur- 
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face density at R, and S cr jt is the critical density for 
lensing which depends on the angular diameter distances 
of the lens and the source. To achieve high signal-to- 
noise S04 stack the lenses. It is thus desirable to re- 
move the redshift dependence of the signal, which comes in 
through S cr it . This dependence is removed if one considers 
AS = 7t x X cr i t . This is the quantity that S04 measure. 

S04 deprojected their AS to obtain the actual 3D galaxy- 
mass correlation function. The deprojection is done via an 
Abel inversion: 

, , x 1 f°° d R f dA ^ „ AS \ 

where r is the 3D separation, R the projected separation, 
and p the mean density of the universe. S04 assumed 
for the present-day matter f^o = 0.27 to obtain the mean 
density p appearing in the inversion. As discussed, all of 
our simulations use fio = 0.3. Thus when comparing to 
S04 3D results we show their results rescaled by a factor 
of 0.27/0.3. 

We calculate AS in simulations by selecting objects 
from the entire halo catalog of a certain simulation in 
accordance with the magnitude distribution of S04 in a 
given band (see their Figure 5). The contribution of each 
object to the correlation function is weighted similarly to 
the weighting of S04. The weight for each lens depends on 
E cr ; t , and thus the redshift of the lens, and on the num- 
ber of background galaxies available for each lens. For a 
flux-limited sample, this weight as a function of redshift 
translates into an effective weight as a function of r-band 
luminosity; to mimic this selection we thus apply the effec- 
tive weight for the S04 sample to our simulated catalogs, 
using a function supplied by E. Sheldon. 

We first compare our results to the direct observable, 
AS, which is commonly referred to as the term (projected) 
correlation function. In the case where the 3D correlation 
function is a power law, AS is the projected correlation 
function modulo a factor which depends on the slope of 
the 3D power law and the average matter density. We 
also calculate the 3D galaxy-mass correlation function. In 
simulations this function can be measured directly rather 
than through inversion of AS by using Eq. (3). 

7.1.1. Luminosity dependence 

We compare the projected and the 3D correlation func- 
tions, AS and £ gm , measured in our simulations to the 
measurements of S04 in the left and right panels of Fig- 
ure 3, respectively. Note that our faintest bin is not exactly 
the same as that of S04, but the faintest galaxies that are 
missing from our sample have only a small contribution to 
the total signal (E. Sheldon, private communication). 

In agreement with S04, we find that the amplitude of the 
correlation function increases with luminosity on interme- 
diate (0.1 — lh^ 1 Mpc) scales, while it is nearly indepen- 
dent of luminosity on larger scales. This is indicative of 
the increase of the effective slope with luminosity. In the 
low luminosity bin our agreement is better for £ gm than 
for AS. As shown by Eq. (3) the value of the 3D correla- 
tion function at given r depends on the behavior of AS 
and of its derivative on all scales above R = r (but with 
contributions closer to R — r more heavily weighted). 

As discussed in §4, we introduce an observationally mo- 
tivated scatter in the luminosity-maximum velocity re- 
lation obtained by matching the luminosity and velocity 




r (hr 1 Mpc) r (h -1 Mpc) 

Fig. 4. — Comparison of the simulation galaxy-mass correla- 
tion function (lines) with the observations of Sheldon ct al. (2004) 
(points) in the u, g, i, and z SDSS bands for the fainter sample. The 
numbers in parentheses denote the actual faint magnitude limit of 
the Sheldon et al. (2004) sample in the cases where it does not co- 
incide with the faintest magnitude for which objects arc resolved in 
the simulations. 



functions. In the right panel of Figure 3 we plot the results 
obtained when no scatter between luminosity and veloc- 
ity is used (thin lines). The net effect of the scatter is 
to reduce the amplitude of the correlation function by a 
factor of <~ 2 — 3 in the high luminosity bin, whereas the 
results for the low luminosity bin are essentially the same, 
cither with or without scatter. This is expected given the 
flatness of the luminosity function at low luminosities, and 
its steep fall at high luminosities. The behavior for AS is 
similar. In what follows, only results obtained with scatter 
are shown. 

The left panel of Figure 3 also shows the difference in 
AS(i?) for two simulations with different values of as- Al- 
though in § 6 we showed that this effect is smaller than 
the cosmic variance for boxes of this volume, the two sim- 
ulations compared here can estimate this effect accurately 
as they used the same random realization of the initial 
conditions and thus arc not subject to differences due to 
cosmic variance. 

We repeated the analysis in the r-band for the other 
SDSS bands and show the corresponding £ gm comparisons 
for the u, g, i and z bands in Figure 4 for the faintest S04 
sample. The luminosities in other bands were assigned 
using the observed density-color correlation for the SDSS 
galaxies (see § 4). The figure shows that the overall agree- 
ment of our simple luminosity scheme, even in the bluer 
bandpasses, is surprisingly good. The agreement for the 
intermediate luminosity bin of each band is similar. 

7.1.2. Dependence on color 

As described in § 4, we have assigned colors to our halos 
based on observational correlation of color and galaxy den- 
sity. In Figure 5, we compare the correlation functions of 
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Fig. 5. — The simulation galaxy-mass correlation function for red 
(g — r > 0.7, dashed line) and blue (g — r < 0.7, solid line) objects. 
The corresponding observational results of Sheldon ct al. (2004) are 
also shown by triangles and squares. For clarity, the results for the 
red objects arc displaced upward by a factor of 10. 



red and blue galaxies separately to those measured by S04. 
In each case, the color separation is made at g — r — 0.7, 
which produces two subsamplcs with roughly comparable 
numbers of objects (as does the same selection in the data). 
The agreement with the data is quite good in both sub- 
samples, especially if we take into account the fact that 
the simulated sample is missing both the dimmest and the 
most luminous objects of the S04 total sample, and that 
cosmic variance is still an effect given the volumes probed 
here (note that we do not plot results from the larger box 
here because it is not as complete at the low luminosity 
end as the simulation shown, but the agreement with the 
larger box is even better). The agreement indicates that 
the simple galaxy density-based color assignment in our 
simulations is sufficient to explain the galaxy-mass corre- 
lations in the data. 

7.1.3. The M{y£.-L scaling 

To quantify the AS-luminosity scaling, many observa- 
tional studies use the "aperture mass" , Mas (Fischer et al. 
2000; McKay et al. 2001; Smith ct al. 2001). This is defined 
as the quantity obtained by two-dimensional integration of 
the A£(i?) up to a certain R. Following previous studies, 
we adopt the outer scale of R = 260h~ 1 kpc (physical). In 
observational studies this scale has been chosen because it 
is small enough to avoid a strong contribution from the lo- 
cal environment of the lenses and large enough to measure 
a significant signal. 

The aperture mass simulation results derived in the r- 
band are shown in Figure 6. The data is divided into five 
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Fig. 6. — The "aperture mass", Mae, as a function of r-band 
luminosity. Solid squares and triangles correspond to the M^e ob- 
tained by fitting the Sheldon ct al. (2004) data with a NFW and a 
SIS profile, respectively. Open squares and triangles correspond to 
results obtained from simulation when fitting these profiles. For a 
more direct comparison to the results obtained from the S04 data, 
the pentagons show the simulation results rc-binncd to match closely 
the S04 luminosity bins. All results are derived by using simula- 
tion ACDMgQb, with the exception of the highest luminosity bin for 
which the ACDM120 run is used to improve the statistics. 



bins with comparable number of objects each. For each 
luminosity bin, we fit the Singular Isothermal Sphere (SIS) 
and the Navarro, Frcnk and White (Navarro et al. 1996, 
1997, hereafter NFW) profiles to the A£(i?) in simulations 
and S04 observations. The selection and weighting of our 
objects is the same as in S04, but the binning differs: S04 
binned their galaxies based on the signal-to-noise ratio, 
whereas our bins are chosen to have an equal number of 
objects per bin. For a more direct comparison we also re- 
bin our objects to match the two lower luminosity bins of 
S04 (pentagons). 

In all cases the NFW and SIS fits give comparable Mas • 
For the SIS fit, AS(i?) = S(i?) (up to an arbitrary con- 
stant), but this is not true for the NFW fit. Note that 
Mas here refers to the integral of AS for both fits. For 
the SIS profile, A/as is the true projected mass (minus a 
mass sheet), whereas the interpretation for the NFW pro- 
file in terms of an actual mass is more complicated, and 
one should perceive Mas simply as the normalization of 
AS rather than a true projected mass, a point emphasized 
by Guzik & Seljak (2001). A simple interpretation of Mas 
as projected mass becomes even more misleading if we take 
into account that for many galaxy-size halos 260/i~ 1 kpc 
is larger than the virial radius. 

Figure 6 shows a good overall agreement in the Mas — 
— L r relation with observations. Differences between the 
aperture masses of simulations and observations are tracked 
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Fig. 7. — The "aperture mass" -to-light ratio, Mas/L t , for the 
entire sample (triangles), for blue (squares) and red (circles) galax- 
ies vs average luminosity of each sample. Simulation results are 
connected by the dashed- dotted line, while the values obtained by 
fitting the AE measured by Sheldon et al. (2004) with the NFW 
profile are connected by the dotted line. Slightly different sample 
definitions arc responsible for offsets in the x-axis. 

back to the bins not being exactly the same in the two 
cases. More specifically, the Mae measured in simulations 
is somewhat larger than observed for the lowest luminosity 
bin. This is due to the difference in the faint limits of our 
and S04 samples: M r — — 18 and M r = — 17, respectively: 
the contribution from galaxies with — 17 < M r < — 18 de- 
creases the Mae and the amplitude of A£(ii) somewhat, 
as previously mentioned. 

The simulations reproduce the Mae oc L r scaling of 
the aperture mass with luminosity found by McKay et al. 

(2001) reasonably well. Interestingly, our results suggest 
that the correlation function flattens significantly for L r < 
10 10 h~ 2 L Q . This is in agreement with McKay et al. 

(2002) , who found that most of the variation of mass with 
luminosity is seen for L r > 1.5 x lO lo /i~ 2 -L galaxies. As is 
discussed below, such a shallowing is also observed for the 
host halo mass-to-light ratio (van den Bosch et al. 2003; 
Vale & Ostriker 2004). 

Using this aperture mass measurement, McKay et al. 
(2001) calculated mass-to-light ratios for the entire sam- 
ple and for samples separated by galaxy type, and found 
that the Maz/L ratio in the red wavebands is insensi- 
tive to morphology. In Figure 7 we compare the Mae / L r 
calculated for our sample with that obtained by fitting a 
NFW profile to the S04 sample (this corresponds to the 
masses designated by solid squares in Figure 6). In addi- 
tion to the value for the entire sample, we also show values 
for the blue and red subsamplcs defined with the color cut 
of g — r = 0.7 (sec § 4 and § 7.1.2). The figure shows that 
the M^/L r ratio is fairly constant and thus insensitive 
to the color of the galaxies. The average luminosities of 
the red and blue simulation subsamples differ from that of 
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Fig. 8. — Top panel: The galaxy— galaxy, the galaxy— mass and 
the mass-mass correlation functions. For the galaxy-galaxy and 
the galaxy— mass correlations galaxies are selected according to the 
SDSS selection function. Only objects with -22 < M r < -18 
are used, close to the flux limited samples of Sheldon et al. (2004) 
and Zehavi et al. (2004). Bottom panel: Inverse cross-bias, bx , 
as estimated by Sheldon et al. (2004) (points) and as measured in 
our simulations (solid line). For clarity, error bars in the simulation 
results are plotted only for scales at which there appears to be some 
discrepancy between simulations and observations. 



the observational subsamples since in simulations we miss 
both the dimmest and the most luminous objects. Note 
that for the data of S04 we find approximately 50% higher 
values of these ratios than were found by McKay et al. 
(2001) for their sample. Our results differ from S04 less 
than the aperture mass-to-light ratio for their blue and red 
subsamples differ. 

7.2. Bias 

The galaxy bias is a measure of how well galaxies trace 
the underlying mass distribution. It is often characterized 
by the bias parameter, b, defined in terms of the galaxy- 
galaxy and the mass-mass correlation functions (Kaiser 
1984) 



(4) 



The correlation coefficient, r, relates the galaxy-mass cor- 
relation function to the mass-mass and galaxy-galaxy cor- 
relation functions (Pen 1998): 

£gm 



V?mm<;ggj 

The mass-mass correlation function is not directly mea- 
surable, but weak lensing observations of the galaxy-mass 
correlations can be combined with measurements of the 
galaxy-galaxy correlation function to estimate the cross- 
bias parameter, 

C gff b 
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Figure 8 compares the inverse cross bias, b~ l , measured 
from simulations with that measured by S04 for a flux- 
limited sample. S04 use the inverse galaxy-mass bias be- 
cause their £ gm is much noisier than £ gg . Their flux- lim- 
ited sample, used to obtain £ gm , spans the r-band magni- 
tude range — 24 < M r < — 17. For £ gg they use the results 
of Zehavi et al. (2004) for a flux-limited sample with r- 
band magnitudes in the range —22.2 < M r < —18.9 and 
selection function similar to that of S04. S04 corrected 
for the difference in average luminosity of the two samples 
by rescaling the correlation length of Zehavi et al. (2004) 
using its scaling with luminosity in 2dF (Norberg et al. 
2001). For the simulation results we use objects in the 
range — 22 < M r < —18 for both £ gm and £ gg and select 
objects to mimic the selection function of S04. The figure 
shows that the cross-bias in simulations is in very good 
agreement with observations. Note that, again, we are 
missing the dimmest and the brightest of the objects used 
by S04. The cross-bias is approximately scale-independent 
with b x s=a 1. 

Figure 9 shows the correlation coefficient r and bias b 
for two volume-limited simulation samples of different lu- 
minosities. In the top panel we also show comparisons 
of the cross-bias, b~ x , for the volume- limited sample of 
S04. The simulation results agree with observations within 
the error bars. A small offset of the simulation results 
toward smaller values of b x may be due to objects with 
—23.0 < M r < —22.5, which are not present in our simu- 
lations due to the small box size and which could enhance 
the clustering signal somewhat. 

Semi-analytic galaxy formation models suggest that b x 
is a direct measure of the bias b on large scales, or, equiv- 
alently, that the correlation coefficient r approaches unity 
for r > Ih- 1 Mpc (Guzik & Seljak 2001; Berlind & Wein- 
berg 2002). Figure 9 shows that the correlation coefficient 
is indeed approximately unity on scales > lh^ 1 Mpc in 
our simulations. This therefore confirms that on these 
scales the cross bias is expected to be a fair measure of 
the standard bias b for both bright and faint samples. 

The scale dependence of bias measured here is in agree- 
ment with previous studies (Colin et al. 1999). Note that 
the lower luminosity sample exhibits an "antibias" at r < 
'ih^ 1 Mpc. For the high-luminosity sample the shape of 
the bias profile is similar to that of the low-luminosity sam- 
ple: it is scale-independent on large scales and decreases 
on smaller scales. The decrease on small scales has been 
attributed to dynamical friction and tidal destruction pro- 
cesses in the high density regions of clusters and groups 
(Kravtsov & Klypin 1999; Zentner et al. 2004). An in- 
crease at small separations, r < O.bh^ 1 Mpc, is likely due 
to the tendency of the brightest galaxies to lie near the 
center of groups and clusters (Seljak 2000; Weinberg et al. 
2004). For the same reason the correlation coefficient of 
the bright subsample increases above unity for separations 
r < l/i^Mpc. 

8. DISCUSSION 

8.1. Luminosity dependence of galaxy-mass correlations 

In the previous section we showed that simulations can 
reproduce the luminosity dependence of the galaxy-mass 
correlation function observed for the S04 sample, if a rea- 
sonable amount of scatter between V max and luminosity is 
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Fig. 9. — Bias for the volume-limited samples. Top panel: Inverse 
cross-bias, b^ 1 , as estimated by Sheldon et al. (2004) (crosses) for a 
sample with —23 < M r < —21.5 and as measured in our simulations 
(solid line) for objects with M r < —21.5. For clarity, error bars in 
the simulation results are plotted only for scales at which there ap- 
pears to be some discrepancy between simulations and observations. 
Middle panel: Correlation coefficient as a function of scale for two 
volume-limited simulation samples with M r < —18 (solid circles) 
and M r < —21.5 (open squares). Bottom panel: Bias as a function 
of scale for the same samples as in the middle panel. 



assumed. The luminosity dependence of these correlations 
is weak at large scales (> l/i -1 Mpc) and strong at small 
scales (0.1 < R < lh^ 1 Mpc). To understand the origin of 
this behavior, we compute contributions to the correlation 
function from host halos ( "central galaxies" ) and from 
halos classified as satellites (see § 3). Figure 10 shows the 
results for the simplest case of the volume- limited samples 
with r-band luminosities in four different bins, [—19, —18], 
[-20, -19], [-21, -20], and [-22, -21]. We adopt the con- 
vention of including the —1 term appearing in the defini- 
tion of the total correlation function in the host halo term. 

The bottom panels of Figure 10 show that the rela- 
tive contribution of the host halos compared to that of 
satellites is considerably larger at high luminosities, as 
expected. The central galaxy contribution at scales r < 
l/i -1 Mpc changes with luminosity only in amplitude, not 
in shape. At the high luminosity end, most galaxies are 
central galaxies and the increase of the correlation ampli- 
tude with luminosity reflects the increase in average mass 
and density of their host halos (Berlind & Weinberg 2002). 
At low luminosities, and at intermediate-large separations, 
the satellite contribution is significant and leads to the de- 
velopment of a 'shoulder' in the overall correlation func- 
tion, since outside the typical satellite radius, one starts to 
measure the density of the environment in which the satel- 
lites are embedded. The existence of this feature makes the 
total correlation function at relatively low luminosities ef- 
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Fig. 10. — Top row of panels: Volume-limited halo-mass corre- 
lation function for four different r-band luminosity bins and con- 
tributions to the correlation function from satellite and host halos. 
Upper left panel: The satellite contribution. Middle upper panel: 
The contribution of host halos. Right upper panel: total correla- 
tion function. Lower left panel: The (total) halo— mass correlation 
function (solid line) for the —19 > M r > —18 luminosity bin, the 
satellite contribution (dotted line), and the host halo contribution 
(dashed line) for the dimmest luminosity bin. Lower right panel: 
The same as in the lower left panel but for the brightest luminosity 
bin. 

fectively shallower than that at high luminosities. 

8.2. Average halo masses and density profiles 
from weak lensing 

An important question about weak lensing observations 
is whether they will allow us to extract information about 
the average halo mass and mass profile of the lenses. Stack- 
ing the lenses to achieve high signal-to-noise makes the 
interpretation of the average mass difficult. It has been 
noted (see, e.g., Guzik & Seljak 2001; Yang et al. 2003) 
that complications arise for several reasons including the 
scatter in sizes and masses of objects and the fact that 
a fraction of galaxies may be located within a group or 
cluster of galaxies. 

These complications could be limited if only the central 
galaxy of each halo was used to measure the correlation 
signal. However, even in this case the interpretation of 
mass is not straightforward. In Figure 11 we show the dis- 
tribution of virial masses 4 for central galaxies in four lumi- 
nosity bins. The distributions in cases with and without 
scatter are shown. The distribution of masses in the case 
without scatter has a finite width because there is scatter 
between virial mass and V max . The figure shows that intro- 
ducing scatter between V^ax and L r significantly broadens 
the mass distributions and makes them non-gaussian. The 
broadening is largest in the most luminous bin. As a re- 

4 Recall that the virial mass is defined as the mass within the radius 
at which an overdensity of 180 times the mean matter density of the 
universe is reached. 
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Fig. 11. — Halo virial mass distributions (defined using an over- 
density of 180 with respect to the mean density of the universe) for 
central galaxies in four luminosity bins with and without scatter in 
the Vmaz-luminosity relation (solid and dotted histograms, respec- 
tively). M^ r s is the virial mass obtained by fitting the AS of the 
corresponding subsample of galaxies, in units of 10 12 /i _1 Mq. For 
the distributions with scatter, the mean and median virial masses 
are indicated with dashed and dot- dashed arrows, respectively. 



suit, it becomes difficult to assign a specific mass to a given 
(even narrow) luminosity range. For example, the figure 
shows that the mean mass is significantly larger than the 
median for every luminosity bin. The best fit virial mass, 
M^ r s , obtained by fitting AS for the corresponding galaxy 
sample to a NFW profile quoted in the figure is a different 
effective averaging, and gives values between the median 
and the mean masses. 

In agreement with Yang et al. (2003), we find a tight 
correlation between M^ r s and the aperture mass Mae- 
Both quantities are derived from the same fit, and thus 
a correlation is expected, but without information on the 
halo concentrations, one cannot be derived from the other. 
Nevertheless, as Yang et al. (2003) showed, this correlation 
is recovered easily if one assumes the well established cor- 
relation between the virial mass and concentration (Bul- 
lock et al. 2001), which in the absence of scatter between 
concentration and virial mass renders the NFW effectively 
a one parameter fit. Given the averaging done to obtain 
AS, the average concentration-mass correlation is suffi- 
cient. In agreement with Yang et al. (2003), we find the 
following scaling for both central galaxies and galaxies in 
low density environments 



M^f = (0.31 ± 0.01) x M AE . 



(7) 



Here both masses are in units of 10 12 ft- _1 Mq. The scaling 
is exactly the same when no scatter is used, since this 
relation reflects the properties of the average mass density 
profiles of our dark matter halos. 

The existence of significant scatter between luminosity 
and mass even for central galaxies indicates that the mass, 
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, derived from weak lensing cannot be interpreted in 
a straightforward way as the mass for galaxies of a given 
luminosity. Instead, there is a rather broad distribution of 
masses. The mass-luminosity relations should therefore 
be quoted and interpreted with caution. Figure 11 shows 
that the amplitude of the relation will depend on which 
mass one considers. The actual distribution depends on 
the still uncertain details of the scatter between mass and 
light, but even if the true relation is different from what 
we have assumed here, it will be incomplete if the dis- 
tribution of masses is not specified. In conclusion, the 
mass distribution at fixed luminosity can be wide enough 
so that the concept of a mean mass becomes obscure. Re- 
ferring to the mass of galaxies that we find is 
w 2 x lO 12 /^ 1 M© as shown for the -21 < M r < -20 bin 
in Figure 11 and is in general agreement with the mass 
found by Yang et al. (2003) and the best fit values derived 
by Guzik & Seljak (2002) from the fit to the SDSS weak 
lensing data. 

A related question is whether a meaningful average den- 
sity profile of halos can be recovered from lensing data. 
Our results show that here, again, one has to focus on 
either the central galaxies or galaxies in low density en- 
vironments. As we show in §8.1, for luminous galaxies 
the galaxy-mass correlation is dominated by the contri- 
bution of central galaxies, especially at small scales. For 
high-luminosity samples it is therefore meaningful to in- 
terpret AE(i?) in terms of the average halo density profile 
at r < 100- 200/i- 1 kpc. 

9. CONCLUSIONS 

We have presented a detailed comparison of galaxy- 
mass correlation functions measured in high-resolution cos- 
mological simulations of the concordance ACDM cosmol- 
ogy with the most recent weak lensing SDSS observations. 
We found that the simple recipe of assigning a luminosity 
to a halo of a certain maximum circular velocity by match- 
ing the subhalo velocity function in simulations to the ob- 
served luminosity function leads to good agreement with 
the observed galaxy-mass correlation and its dependence 
on luminosity, if an observationally-motivated amount of 
scatter is introduced in the V max — L r relation. Our main 
results and conclusions can be summarized as follows: 

• The simulations reproduce the galaxy-mass corre- 
lation function measured by Sheldon ct al. (2004) in 
SDSS and the observed dependence of its shape and 
amplitude on luminosity. Interestingly, the agree- 
ment for bright samples seems to require some scat- 
ter between luminosity and circular velocity. The 
amount of scatter required is consistent with the 
observed scatter. 

• The galaxy-mass correlation function for central 
galaxies is a strong function of the galaxy luminos- 
ity (halo velocity), while £ gm for satellite galaxies is 
only weakly sensitive to luminosity. The luminosity 
dependence of the correlation function as a whole 
is thus determined primarily by the increasing con- 
tribution of bright central galaxies relative to the 
satellite galaxies at bright luminosities. Conversely, 
the correlation function gets shallower with decreas- 
ing luminosity because it is increasingly dominated 



by the contribution from the satellite (non-central) 
galaxies in halos at intermediate and large scales. 

• We use the color-density correlation observed in the 
SDSS to assign colors to the galactic halos in sim- 
ulations and compare the simulation results to ob- 
servations in other SDSS bands, from u to z. In 
each band the agreement between galaxy-mass cor- 
relation function in simulations and observations is 
remarkably good. The simulations also reproduce 
the observed trend of the galaxy-mass correlation 
function with the g — r color. 

• In agreement with previous studies, we find that 
the aperture mass-to-light ratio, Mas/L v , is in- 
dependent of galaxy color. For L r > 10 10 h^ 2 L Q , 
M&y, I L r is approximately independent of luminos- 
ity. For both central galaxies and galaxies in low 
density environments the best fit virial mass M^J 1 
correlates tightly with the observed aperture mass. 
The best fit virial mass lies in between the median 
and the mean actual virial mass when scatter be- 
tween Vmax and luminosity is used. 

• We compare the cross bias, b x = b/r, measured in 
simulations and observations and find a good agree- 
ment at all probed scales. The correlation coeffi- 
cient, r, obtained for a volume-limited sample sim- 
ilar to that analyzed by S04, is close to unity on 
scales > Ih^ 1 Mpc. This indicates that the cross 
bias measured in weak lensing observations does 
measure the actual bias b of galaxy clustering on 
these scales. 

• We show that for luminous galaxies (M r < —21) 
the galaxy-mass correlation function at r < 100 — 
200h~ 1 kpc can be interpreted as the average den- 
sity profile of these galaxies. We also show that 
the masses obtained by the AS(i?) measured in 
weak lensing observations cannot be interpreted in 
a straightforward way as the mass for galaxies of a 
given luminosity. In the presence of scatter between 
mass and light, the obtained mass- luminosity rela- 
tions should be quoted and interpreted with cau- 
tion. 
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